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Abstract 

A representation is describee whose basis levels combine the important phy- 
sical aspects of a finite set of plane waves with those of a set of Bloch tight- 
binding levels. The chosen combination has a particularly simple dependence on 
the wave vector k within the Brillouin Zone , and its use in reducing the standard 
one-electron band structure problem to the u;iual secular equation has the advan- 
tage that the lattice sums involved in the calculation of the matrix elements are 
actually Independent of k. For systems with complicated crystal structures, for 
which the Korringa-Kohn-Rostoker (KKR) , Augmented-Plane Wave (APW) and Ortho- 
gonal! zed -Plane Wave (OPIV) methods are difficult to apply, the present method 

leads to results with satisfactory accuracy and convergence. 
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1 . Introduction 

The method described below evolved from an attempt to obtain the band struc- 
ture of a system such as molecular hydrogen in a relatively complex crystal struc- 
ture, and in a variety of density regimes. For certain regions of the density it 
is expected on general grounds that neither the low density tight-binding approach 
(with a representation of linear combinations of atomic orbitals (LCAO) Blocb 
functions) , nor the methods using a representation with a basis of simple plane 
waves (PW) , are physically adequate. 

For reasons principally connected with the structure, the other standard 

1,2 

methods are also not entirely adequate. ’ The Korringa-Kohn-Rostoker (KKR) and 
Augmented-Plane Wave (AP^V) methods not only require an extraordinary amount of 
computational effort, but are based on a muffin-tin approximation to the actual 
one-electron potential. * ’ _ This mea ns a "sphericalizatlon" (taking the average 
over angles) of the potential arising from the contents of a unit cell, a pro- 
cedure which is difficult to justify when the molecules in the crystal have no 
obvious spherical symmetry. Although such models yield useful physical informa- 
tion especially at lower densities, it is difficult to estimate their accuracy, 
particularly at higher densities, where steric effects and the requirements of pro- 
per crystal symmetry may become important. The effects of the latter on the resulting 
band structure may well be important. 

Furthermore, there is often no clear cut separation between core levels (actually 

non-existent for hydrogen) for which tight-binding is adequate, and the rest of 

band levels (valence and conduction) , which would make an orthogonalized plane 

wave (OPW) method f^p^oprlate. Even if one makes an arbitrary separation between 

valence and conduction levels, and treats the first with tight-binding and the 

0 

second with OIAV functions, orthogonalized to the valence levels, one still has 
the possibility of significant overlap of these "core" levels in situations where 
large variations in density are of physical ii>terest. 
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For these reasons It Is natural to Investigate alternative representations 

whose basis levels combine in some way the advantages of both the LCAO functions 

(with their physically correct atomic behaviour near the nuclei) and the PW which 

are more satisfactory in the region between atoms. One such basis set was recently 
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used by Ramaker et £l in exact-exchange crystal Hartree-Fock-calculations of mole- 
cular and metallic hydrogen. Another, based on a more general and flexible approach 
is described below. It is a modification of an idea used successfully by Brown 

g 

and Krumhansl, which was shown to be mathematically equivalent to but in fact 
more general than the orthogonalized-plane wave method. 

In the next Section the representation will be developed and its basic pro- 
perties described. Section III is devoted to a discussion of the application 
of the representation to the solution of the one electron problem in crystals. 

In -Sec t ion - IV w e di scttss-fartefly- t h e po ssible applications of the method, which 

although originating from the study of a specific material is of more general 
interest and can be used to study the electronic structure of a wide class of 

materials. The specific results for the band structure of molecular hydrogen will 

9 

be given in a following paper. 

II. Representation 

The representation we introduce is formally incomplete but only in the 
sense that it has a finite set of basis wave functions. This set is made up of a 
finite number of PW and a set of specially constructed Bloch levels. It is con- 
structed in such a way that the whole set is orthonormal, and although the set 
finite, linear combinations of them are expected to give variationally good approxi- 
mations to the eigenfunctions and corresponding eigenvalues. This expectation is 
based on the physical way the set is constructed, which will be explained in what 
follov s. 

Consider first a monatomic (for example, a simple cubic) lattice with lattice 
constant n and LCAO-Bloch level h (r) defined with atomic orbital $(r) as 

*— K ^ 


follows: 
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h (r) = — ) 

^ /N R 


( 1 ) 


where N Is the number of cells in a volume fl, R designates their position vectors, 
and k is the Bloch wave vector. Expressing this Bloch level in its well known form 




•I 


yn-K 


k-K 
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where K is the set of reciprocal lattice vectors corresponding to R , it is 
easy to see that 



where $ is the Fourier transform of $(r). 

a, 

For the purposes of defining a trial function, may be any localized 

orbital, and not necessarily an atomic one. This observation will be used to con- 
struct a particularly convenient type of Bloch level. But instead of defining 

it directly (l.e. in r-space) it is inferred from conditions imposed on c . In 

q 

A 

this way it is easier to enforce( through them) the properties that one would like 
the Bloch levels to have. First, some general observations: 

One expects the eigenfunctions not to change too much very near (and parti- 
cularly Inside, if there is a core) the atoms or molecules forming the solid from 
the values they assume in corresponding free atoms or molecules. This remains true 
even at fairly high densities. Thus one wants to include in the basis set, Bloch 
levels built with at<imic or molecular orbitals to obtain a good representation in this 

region. But it is clear that for this purpose only those components c. with 

K"K 

sufficiently large K are relevant (here k is assumed to be restricted to the first 

Brillouin zone B ). On the other hand, if the itinerant or free electron charac- 
o 

ter becomes important (as it u'll at high densities) plane waves with wave vectors 
(about the origii^ not too large in terras of 2n/a are obviously Indicated. We now 
construct basis functions incorporating these features. The Bloch level is first 
modified by truncating its Fourier components of low wave vectors, say G in some 
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flnlto subset G of the reciprocal lattice K. In tnis way the plane waves with 
wave vectors k-G have been set free to be included in the basis net as inuepondent 
members orthogonal to the Bloch levels. (For simplicity in some of the algebraic 
manipulations the subset G may be chosen symmetrically to include both G and »G, 
although this is not required in general by the method.) For the simple cubic lattice 
case, for example, we may choose G to be the set of all reciprocal lattice vectors 
within or on the surface of a cube centered at the origin, and with faces perpen- 
dicular to the axes. Further, let T be the complement of G, that is GoT is empty 
and GoT = K. Next, the Bloch levels of the basis are to be chosen to have as simple 
a form as possible, a requirement for both analytical and computational purposes, 
in particular, the most simple functional dependence on k is essential. 

In the case of a Bravais lattice, a set of Bloch levels satisfying these 
criteria can be taken to have components 

Y TO L h (rs \ <^> 

K -2 “ ~ 


where the characteristic function given by 




1 , if X e A , 

1^0 , otherwise. 

Here $(r) is a localized orbital. Figure 1 shows a schematic one-dimensional 
example of the procedure just outlined; there the dotted curve represents the Fourier 
transforms $ of a localized orbital and the dipontinuous curve the components 

q 

(n/N)^c^ given by Eq. (4); note also that the set G contains by choice only the 
reciprocal lattice vectors 0 and i 2n/a. 

The levels defined by Eq. (4) all have the properties of Bloch levels, and 
can of course be written as 


h.(r) = “ y C 

K, 


i(q-K) .r 

0 


< 5 ) 
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Thls reduces, for JJ, = ® to the standard form; 



ik.r 


0 ^ 



IK.r 

e 


(6a) 


and is equivalent also to: 




ik.r 


/T 


e 


- R G6G 


(6b) 


where the quantity in [ ) clearly has the periodicity of the lattice. The pre- 
factor in the expression for c is not Important except to keep track formally, 

%> 

and in a consistent way, of the various constants and factors involved. (It 
cancels when normalizing the functions.) 

The norm of ^^^^(r), Hht] , is Independent of k and is given by 



(7) 


or equivalently by 

^ 0(r)|$(r-R)> “ 7 • 

„ GeG 

K /w 


with the normalized levels (h (r^)/|lh||) the corresponding Wannler 


function w(r) can be obtained, and is given by 


w(r) 


=2 __i V c 

llhliyiffi L ft 

a. 


(9) 


which in this form is automatically normalized. It is, of course, orthogonal 
to w(r-R) for 1^0. Substituting in Eq. (10) for c , one gets the more explicit 

^ Q 

% 

form: 
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or 


w 


' h 


V N ^ 

5 - hv-r) - 1 


i 


ft e 
G 




w ( , (lOb) 


where for the case of a simple cubic Bravais lattice 

1 


0 


ik*r 



sin(nx/ja) sin(TTy/a) 

(TTx/a) (TTy/a) 


sin(TTz/a) 

(nz/a) 


( 11 ) 


is the empty lattice lowest band Wannier function. 

It is clear from the form of h, (r) and w(r) that these functions have the 

k ~ ~ 

right behaviour near and at the lattice sites R, particularly if the finite set G 
does not contain large wave vectors. And for all jG c G, h^(^) Is automatically ortho- 

r*i/ 

gonal to the plane waves with wave vector Is-jG. 

In this way we have an Incomplete but L»rthonormal basis set which v/ould clearly 
■.be sufficient for a monatomic lattice if it were not necessary to use more than one 
localized $(£,)• 


Except for small k, the Bloch level hj^(jr) just defined will in general not be a 
good approximation to the solution Y|^(_r) of the one-electron problem of the crystal, 
if G is empty (l.e. if no PW are Included in the basis). The functions 
h^(r) and can differ substantially for larger k, particularly near the bound- 

arles of the Brlllouin zone, simply because the Fourier components of exp(-ik.r) \Jt (y) 
are functions of k, while those of exp(-lk.r) h, (r) are not. Nevertheless, 
considering their expansions in reciprocal space, as K increases the differences in 


their components decrease, since by construction both functions have the same form 
inside the atoms. Therefore, by truncating the components of low K, and including 
the corresponding PW with wave-vector k-K in the basis, we will increasingly improve 
the approximation as the number of PW increases. 
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Certainly It would bo a better approximation to start by truncating the 

TI) 

usual tight-binding Uloch level h|^ (jr) (defined with $(jr))and choosing components 


a 



( 12 ) 


BO that 



(r) 


Q 


I 

KCT 


\-K» 




(13) 


But this would not have the immense computational advantages of form (6), which 
permits all the terras there co be expressed in lattice sums independent of k. 
Nevertheless, for some cases higher accuracy requirements together with the 
necessity to keep the number of PW within reasonable limits might make it mandatory 
to use better Bloch levels than those defined by Eqs. (6). [One way of defining 
these that would still give lattice sums independent of k, is to take 



up to some order, but of course the higher the order chosen the more cumbersome and 

time consuming become the computations.] 

For the case where a set of more than one linearly independent localized 

orbital must be used, a special Bloch level h. (r) must be included for each. If 

k ~ 

the cell con tains several atoms, say M atoms, with position vectors B. (i=l,2, , . . ,M) , 
a set hj^(^-B^) (i=l,2, . . . ,M) of linearly independent Bloch levels, or M indepen- 
dent linear combinations of them, must be included in the basis set. All the special 
Bloch levels are assumed constructed with a truncated set of plane waves of wave 
vectors ^-G with reciprocal lattice vectors G belonging to one and the same subset 
G. The basis will then contain for the same k (other than the truncated set of 

— rw 

plane waves) a set of linearly independent Bloch levels orthogonal to them but not 
in general to each other. An orthogonallzation procedure must then be used to get 
an orthonormal basis set. The use of this orthonormal basis ultimately results in 
a secular equation with the energy eigenvalues residing only on the main diagonal. 
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rid has distinct analytical and computational advantages. The selection of one 
particular linearly independent set of Dloch levels (over other possible equivalent 
sets) depends on a Judicious evaluation (as far as this is possible) of how well 
they represent the true eigenfunctions of the cryf al, and how their form may help 
the orthogonalization procedure in efficiently producing a physically convenient 
orthogonal set. 

Let the initial set of Bloch levels, before the orthogonalization procedure , 

be a set of linearly independent combinations defined by 

M 

^nk^^ ■ Y. ^nj **Jk^~^ » n=l,2,...,M, (15) 

~ J=1 

where the constants a . are df^termlned In a moment. Here the h., (r) are the 

nj ik ^ 

fS0 

Bloch levels defined for simplicity (but without loss of generality) with only one 
localized orbital in one of the monatomic sublattiices of the basis. Hence: 

Sk's’ = •Sc<£.-Sa> • ( 16 > 


Now we use the Gram-Schmldt orthogonalization procedure to get from 

an orthonogonal set {g ) The g . have the following recursion relations 

‘ nk ^ 




Ik' 


~ f. 


le„J! II f 


rak" " nk' 


(IV) 


l®nk>=^„ -,4 - . for n=2,3...,M. 


and the norms II 6^1! are given by 


’nk" 


= 1 




-I 


»=0 "‘nk"' 


( 18 ) 


Tliose may be used in slightly modified form which subsequently reduce 
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tho numerical work. Lot bo exproa,3od first as linear combinations of 

“jk <r> • 


Then 


^‘^nk> = y ' for n=2,3,...,M. (19) 

J 


^mk^^nk^ 

/V 


y y **mik ^nj ^^’ik^^jk^ 
1 J 


( 20 ) 



(Note that in general these are functions of k.) 


Further 


V * 


^nk" 


= U .. 

1 d 


hi ®nj 

(o) 


(22) 


Next let an orthonormal (incomplete) basis set (r) , a SA, k e 

/W — 

be defined by 


(o) 

Wi)' < 


(o) - i(k-G).r for Of= G e G , 

V — * e — 

Si ^ 


(O) 


for a=n,l<n<M. 


Then A = G {n,l<n^} . The superscript zero indicates this is a basis in 
which to expand the unknown variational approximations to the eigenfunctions 


Y^(r), i.e. 


^ (o) 

\'i.> = 2- *ak • 


asA 


(23) 


( 24 ) 
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Equntion (24), p,s an expansion of the ono-cloctron function, will be used in the 
next section as a trial function for the one-electron problem in crystals. Note 
that, although incomplete, the finite basis sot (23) is orthonorraal and contains by 
construction localized orbitals appropriate for the cores of the 
molecules forming the crystal and plane waves adequate for the Intermolccular 
region. Therefore, we can expect linear combinations of them to be good approxi- 
mations for the eigenfunctions of the lower bands, the accuracy improving as the 
number of PW in G increase it ^particularly for ^ near the boundaries of the Brilloulr. Zone. 


Ill , Application to the solution of the one-electron problem in crystals 

Substituting Eq. (24) into the one-particle Schrodinger equation for the 
crystal, the band structure problem reduces to 

I "oBk ' \ *ak all a e A , „5, 

BCA ~ ~ ~ ~ 

with 



(26) 


Here H is the single particle crystal Hamiltonian. The reason why only one k is 
involved everywhere is that H is a linear operator invariant under the translation 
group of the crystal, for which; 



The matrix elements H g. are given by: 

2 

**G’Gk “ 2ra 


G'G 


G’-G 


(28) 


H 


Gnk 


\ *’njk I ^ !l®nk II 


(29) 



- 12 - 


^*n'nk ^n'jk 

X j ^ ^ ^ ^ ^ ^ 


where, as usual, the plane wave matrix clement of the local one-electron 
crystal potential is given by 

“k = n \ • 


(31) 


with 


» j -iK.r 

V„ = dr e ^ V(r) 

K J ~ ' 


(32) 


and 




V(r-R) 


(33) 


Because of the special form ( 6 ) of h., (r) the products <Ch^, Ih^. ^ and the 

ik ~ '' ik' jk'^ 

/ f \ ^ ~ ^ 

matrix elements |Hjh > and ' can be expressed in terms 

of reciprocal (or reciprocal and direct) lattice sums which are Independent of the 
point in the Brillouln zone (all the k-dependence being factored out). For the 
case of only one localized orbital but with a basis of several atoms, we have 


/u lu \ N lk.(B.-B,) 

n ® Sij , 




(34) 


(35) 


and 


<"ikl = n It 


ij 


(36) 


where 


o Y I* l2 1K.(B.-B.) 

®1J = Z. I«K I e ~ ^ ~J 

K0T 


(37) 
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(38) 


id 




|5„ 1" e‘£-<Sl-Sj> 

K 


(39) 


« . » „ iK.B. 

“&i Vii ' ~ 


(40) 


and 


s? 


” = I y 

KfeT kHt 


S'" * tJ -i(K* .B.-K.B,) 

\ \’-K e ^ ~ . 


(41) 


These lattice sums can be expressed in part as direct lattice sums, using the 
convolution theorem or by application of Eq. (6b). For example; 

®1J = I <«<£> 1«<J?-Si-Sh-S» - ttI! <«> 


Gee 


Prom this, S . . and S". . can be obtained respectively by taking tho gradient and 
M.J ij 

the negative of the Laplacian with respect to the spatial variable. A similar 

result can be obtained with S^.and S^. , but here It would be of no advantage if only 

SP ij 

the Fourier transform of the potential Is available. 

The number of different lattice sinns that must be actually computed is greatly 
reduced by exploiting crystal symmetries. First of all, the sums are invariant under 
a transposition of indices, except for (which only elmges sign) and . In 

general a simultaneous change of B . , B. and G (in the case of S ) under the same 
cubic or other symmetry will also leave S^j, ®ij 

take S ' . into the corresponding symmetric vector. In this way, for example, the sixty 




(10) 


four S. . sums of the Pa3 (or a-N_) crystal stmeture are reduced to only four, and 

Xj 2 

s_. sums to only two for each G, and in both classes of sums this leads to an 

gP ~ 
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enomous reduction in computational time. 

Once the lattice sums are evaluated, wo can proceed to solve the secular 
eigenvalue problem (25) tor a particular k by first obtaining the corresponding 
basis set (23) with the help of Eqs. (19)-(22), then the matrix elements H g. 

with Eqs. (28)-(30) , and finally diagonall eing Eq. (25). In this way we obtain 
the valence and lowest conduction bands and the coefficients In the expansion 
of the corresponding eigenfunctions in terms of the basis set (23) . 

IV. Summary and Conclusions 

We have shown that subspace spanned by the orthonormal finite basis set of 
functions (23) , can be expected to be u satisfactory approximation to the actual 
eigenfunctions. The set is of manageable size and at the same time leads to good 
convergence because by construction it contains orbitals which represent well 
both the core and Intermolecular features. This is accomplished in a rather simple 
way, with a few plane waves and orbitals depending on k only through a factor 
exp (ik.^. It leads to lattice sums independent of ^ when calculating the 
matrix elements of the secular eigenvalue problem (25) , to which the band structure 
calculation has ''een reduceo. It is then necessary to evaluate these lattice sums 
only once for a given lattice parameter and crystal structure. Finally, it is 

straightforward to obtain the necessary matrix element for any k even for 
complicated crystal structures. 

The mei.hod is potentially more flexible than the KKR and APW methods, since 
it does not require the muffin-tin approximation for the crystal potential, and 
is readily adaptable to molecular and complex crystal structures. The level of 
analytic complexity and computation difficulty does not exceed that of the 
standaixl methods. When, compared specifically with the OPW method its main 
advantage is the simpler and more flexible formulation, which makes no particular 
reference to core levels. 

There is another sense in which we can also expect good convergence, namely 
the evaluation of the lattice sums (37)-(41). Notice that with the present method 
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the choice between the purely reciprocal space lattice sums or sums divided 
between reciprocal and direct space will depend on the convergence properties of 
the parti'" alar case at hand. For example, in the case of solid molecular hydrogen, 
$(r) can be taken as a Is-orbital, namely: 


$(r) 



.-ar 


(43) 


with Fourier transform: 


a 


/ srra 

I * 


(44) 


For the direct lattice sura in Eq, (42), we then find: 
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<^$(r-r’) l§(r')> = e”®*" (1 + ar + i aV) , (45) 


which leads to rapid convergence in direct space for the S^j and Since 

s and involve $(r) (wh*ch falls exponentally with r) and V(r) (which goes more 
Gi ij ~ ^ 

or less as r ) a similar conclusion can be drawn about their convergence in 

direct space. But note that in reciprocal space, the lattice sums of Eqs. (37)- 

I -4 -2 

(41) also converge faitly well, since 5^^ goes as K , and Uj^ approximately as K . 

rv 

The method has been applied to solid molecular hydrogen in the a-N^ structure 

for densities varying from the zero pressure value up to more than ten times that 

9 

value. The details and results will be published in a later paper. it is 
sufficient here to say that we obtain results with satisfactory accuracy and con- 
vergence. They corroborate well with the results of other methods. However, as 
mentioned earlier the method described above is of more general interest in the 


context of band theory. 
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Figure 1. Schematic one-dlmcnslonal example of components (0/N)^c of a member 
of the new representation given by Eq. (4) (discontinuous curve) In terms of the 


Fourier transform of a localized orbital (dotted curve). The reciprocal lattice 
vectors correspond hero to q/(2TT/a) = integer. Note that c is identically zero In 

” q 

the central zones (corresponding to a choice here of a set of reciprocal lattice 


vectors G = { ~2TT/af 0 , 2TT/a } ) and constant within each zone corresponding 

to the reciprocal lattice vectors falling outside G (set T) . 
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q/( 27 T/a) 





